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In the simple model of a massive vector field in flat space-time, we derive a kinetic equation of non- 
Markovian type, describing the vacuum pair creation under the action of external fields of different 
nature. We use for this aim the non-perturbative methods of kinetic theory in combination with a 
new element when the transition of the instantaneous quasiparticle representation is realized within 
the oscillator (holomorphic) representation. We study in detail the process of vacuum creation of 
vector bosons generated by a time-dependent boson mass in accordance with a conformal-invariant 

■^j- ' scalar-tensor gravitational theory and its cosmological application. It is indicated that the choice 

of the equation of state (EoS) of the Universe allows to obtain a number density of the vector 
bosons that is sufficient to explain the observed number density of photons in the cosmic microwave 

£SJ ' background radiation. It is shown that the vector boson gas created from the vacuum is in a strong 

' non-equilibrium state and corresponds to a cold dust-like EoS. 
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I. INTRODUCTION 



The vacuum creation of massive vector bosons in intense fields of different nature is widely discussed in the literature 
and it is particularly interesting since the massive vector field is the simplest example of a quantum field theory 
(T) | with higher spin and can be used to model aspects of gauge field theories. To mention examples: the massive 
OO ■ vector bosons of the standard model play an important role in different physical problems in cosmology (see, e.g., 
the massive vector field can be considered as a simplified version of QCD 10], and the p meson as a "massive 
photon" plays a key role in the diagnostics of hot and dense nuclear matter in heavy-ion collisions via its decay into 
lepton pairs [Tl|. 

Different methods have been used for the investigation of the vacuum quantum effects for vector particles, e.g., the 
imaginary time method [f| , the Bogoliubov canonical transformation method |(j , and the oscillator representation 
0. The most attention has been devoted to the massive vector particles with gyromagnetic ratio g — 2. The 
probability of pair creation of a vector field with arbitrary g in the constant electromagnetic field has been considered 
O i| in 0. The generalization to arbitrary time dependent electric field has been performed in using the method of 
diagonalization of the Hamiltonian, where problems have been encountered in treating pair production by electric 
fields for the case g — 2. This was unexpected in view of a successful evaluation of the Lagrange function for a 
constant field in the one- loop approximation Q. 

In the present work we give a kinetic description of the vacuum creation of charged massive vector bosons under the 
influence of a time dependent spatially uniform electric field of arbitrary polarization. We also consider the possibility 
■ °f a time dependent mass which represents a new independent mechanism of vacuum particle production. 

The use of kinetic methods allows one to obtain a rather general solution of the non-perturbative problem for an 
arbitrary time dependence of the strong external fields. The non-perturbative approach is particularly appropriate 
for fastly changing fields such as, for example, in the case of time dependent vector boson masses in the vicinity of 
the cosmological singularity, see Sect. II VI 

The construction of a kinetic theory of vacuum particle creation on a dynamical basis requires the time dependent 
quasi-particle representation (QPR) 0,0. We use the oscillator representation for this aim 0] as the most effective 
instrument for the derivation of dynamical equations in the QPR in Sect|n] We introduce here two types of QPR: 
the complete one (based on the full diagonalization of all physical quantities in the Fock and spin spaces) and the 
incomplete one which lets the spin projection uncertain. Further in Sect. IIIII we use these results for the derivation 
of the kinetic equations (KE). A new feature of the obtained system of KE is the presence of a tensor distribution 
function in a rotating coordinate system with the orientation defined by a time dependent kinematic momentum that 
results in a new type of non-Markovian processes. A significant simplification is achieved when the non-Markovian 
effects are neglected. The case of the absence of an electric field is considered in detail when the vacuum creation is 
caused entirely by the time dependence of the mass. The system of KE splits into separate equations for the transverse 
and the longitudinal components which can be investigated numerically. As an application we reinvestigate in Sect. 
IIVI the creation of massive vector bosons in the early Universe within a conformal invariant scalar-tensor theory of 
gravitation as suggested earlier by Pervushin and collaborators 0, 0] . In this approach the time dependence of the 
scalar field entails a cosmological evolution of all particle masses which, according to Hoyle and Narlikar may serve 
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as an explanation for the cosmological red shift alternative to the Hubble expansion. In the present approach, we 
are able to remove the singularity in the density of the produced longitudinal vector bosons which has been reported 
previously • We present a solution of the KE for a toy model where the time dependence of the scalar field is given 
and show that the density of vector bosons created in the early Universe corresponds to the number density of cosmic 
microwave background (CMB) photons. As it will be shown in Sect. II VI the massive vector boson-anti-boson gas 
created from the vacuum is a cold one. It leads to a coarse grained pressure that proves to be approximately zero, see 
Sect. (0. At the same time, the energy density grows at large time scales proportional to the particle mass. That 
leads to a dust-like equation of state (EOS) of massive vector boson - anti-boson gas, in distinction from the static 
stiff EOS of the usual gas of massive vector bosons [13. Finally in Sect. I VII we summarize and present the conclusion. 
We use the metric g^ v = diag(l, —1, —1, —1) and natural units h = c = 1. 



II. THE QUASI-PARTICLE REPRESENTATION 

We consider the vacuum creation of charged massive vector bosons in the flat Minkowski space-time by the action of 
two mechanisms: (i) a time variation of the boson mass m(t) and (ii) the action of some classical spatially homogeneous 
time-dependent electric field with 4-potential (in the Hamilton gauge) 

A^t) = (0,A 1 (t),A a (t),A 3 (jt)), (1) 

where the corresponding field strength is E = —A, and the overdot denotes the time derivative. 

Thus, the field can be considered either as an external field, or as a result of the mean field approximation, based 
on the substitution of the quantized electric field A k (t) with its mean value ( A k (t) ) = A k (t) , where the symbol (...) 
denotes some averaging operation. The time dependence of the vector boson mass can be interpreted as a result of 
the coupling to an average Higgs field. In the kinetic theory, the consideration of fluctuations leads to the collision 
integrals |l6(. Thus, the mean field approximation corresponds to the neglect of dissipative effects. 

We will restrict ourselves to the simplest version of the theory with the Lagrangian 

£(x) = -D^ulD^v" + m 1 u* v u v , (2) 

where _D M = + ieA^, and e is the charge of the vector field including its sign. Eq. J5J) leads to the equation of 
motion 

(£>„£>" + m 2 K = (3) 

with the additional constraint 

= 0. (4) 

The transition to the QPR can be realized in different ways, e.g. by means of the time-dependent Bogoliubov 
transformation Q, or with the help of the holomorphic (oscillator) representation (OR) We choose the OR 

being the simpler method. The OR can be introduced in the spatially homogeneous case and it is based on the 
replacement of the canonical momentum by the kinematic one p — > P = p — eA in the dispersion law of the free 
particle u>(p, t) = \J m 2 {t) + P 2 to be used in the standard decomposition of the free field operators and momenta in 
the discrete momentum space |l7j 

^^ X ) = -'Y J \We- f "{»W(p,t)-l(-)(-p ) t)} 1 (5) 



where V = L 3 and pi = (2tt/ 'L)n,i with an integer rij for each i = 1,2,3. The substitution of the field operators (J5J 
into the Hamiltonian 

H ■ = - J dx (n^ + D*u;Dw^ + m 2 u;^) (6) 
brings it at once to a diagonal form in the Fock space which corresponds to the QPR 



fT = -£>(p,i) 



(+> (p, i)a ( -^(p, t) + (-p, *)&<+>"(-?, t) 



(7) 
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However, this quadratic form is not positively denned. In order to exclude the [i = component with the help of the 
additional condition Q , it is necessary to derive the equations for the amplitudes , . 
Substituting the Eqs. JHJ) into the Hamiltonian equations 



SH 



SH 2 , 
7Th = — - — = m u u 
8u^ 



(8) 



we find the Heisenberg-type equation of motion for the time-dependent creation and annihilation amplitudes 

^(-P,*) = ^A(p,t) a W( Pl *)±iw(p > *)6^(-P.*) J 



where 



A(p,t) 



w(p,t)' 



(9) 



(10) 



Analogous equations were obtained in the work [ll| for the case of scalar QED on the basis of the principle of least 
action. 

Thus, the Hamiltonian formalism in the OR leads to the exact equations of motion JIJJl for the creation and 
annihilation operators of quasi-particles, depending on the "natural" representation of the quasi-particle energy uj(p, t) 
in the external field QJ. 

The additional conditions (@J may be transformed now with the help of Eqs. © to the following form (i — 1, 2, 3) 



^(p,i)4 ±) ( P ,i) = p^ (±) ( P ,i), uj( P) t)b^(-p,t) = -Pi&f^-p,*). 

These equations allow to exclude the fi — component in the Hamiltonian (Q) with the result 

H = "fo *) {«! +) (P. *k H (P. *) + t ] (-P, ^! +) ("P, *) - 



(11) 



1 



c 2 (p,t) 



( P ,t)) (p fc 4- } ( P ,t)) + (Pi^c-p,*)) (p fe &i +) (- P ,t))]l. (12) 



da 



The next step is the additional diagonalization of the quadratic form (|12|) by means of the linear transformations 



a (±) (p,t) = Pa (±) (p,t) = e 1 a^ ±) (p,t)+e 2 a 2 ±) (p,t)+e 3 -4 ±) (p,i), 



b (±) (-P, t) = E(3^>(-p, t) = ei pr> (-P, *) + e 2 /3r j (-P, *) + e3-/3f ; (-p, t), 



?(±) 



(±), 



(13) 



where [ei(p, i), e2(p, t), es(p, i)] determine the local rotating basis built on the vector e 3 = P/|P|. These real unit 
vectors form the triad, 



(14) 



The presence of the factor u/m in the non-unitary matrix E in Eq. Q13|) . leads to a violation of the unitary equivalence 
between the (a, 6) and (a, /?) representations. 

The transformation (|13fl leads to the positively defined Hamiltonian 



P = ^>(p,i) 



a\ +) (p, t)«f > (p, t) + /?r J (-P, tW (-p, t) 



„(-) 



(-), 



?(+)/ 



(15) 



Let us write the equations of motion for these new amplitudes as the result of a combination of Eqs. © and l|13|) 



of (P. *) = ^ A(p, i)/3 4 (T) (-p, i) ± iw(p, t)af ] (p, i) + 77y (p, t)a 3 (p, i), 



-P,*) 



i A(p, t)4 T) (p, t) ± iw(p, i)/?f ) (-p, t) + m (p, i)# (-p, t). 



(16) 
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The spin rotation matrix rjij is defined as 



e!e 2 ^e!e 3 

-eie 2 ^e 2 e 3 
-%eie 3 -77re 2 e 3 -A ro , 



(17) 



where A m = —m/m+ A. Together with the Hamiltonian Q, the operators of total momentum and charge take also 
a diagonal form. However, the spin operator 



Si = E ijk J dx [u* k TT* + TTjUk ~ U*TT* k - 7Tfcl 

has a non-diagonal form in spin space in terms of the operators a^' and 

Sk = itijk ^ 



a[ +) (p, t)a ( r ] (P, t) - (-p, t)$ +) (-p, t) 



(18) 



(19) 



In particular, the spin projection to the momentum 3-axis is 



5*:, 



E 



a[ +) (p,t)a y 2 '(p,t)-a^>(p,t)a\ >(p,t) 



.,(-) 



(+>, 



.,(-), 



+ p { 2 \- p ,t)p[ + \-p,t)-0{ \-p,t)4 + \-p,t) 



(20) 



Thus, this representation can be called an incomplete quasi-particle one with non- fixed spin projection. The operator 
(|20|l can be diagonalized with a linear transformation to the circular polarized waves basis |l7| 

cf\p,t) = R^a^(p,t), 



d^(-p,t) = i^*/?£ ±3 (-p,t), 



with the unitary matrix 



R 



1 Ti 
±i 1 
\/2 



(21) 



(22) 



As a result, the new amplitudes c^', d^' in the QPR correspond to the creation and annihilation operators of charged 
vector quasiparticles with the total energy, 3-momentum, charge and spin projection on the chosen direction 



H(t) = £>(p,t) 
p 

n(t) - E p 
p 

Q = e E 
Ss(t) = J2 



Cj (+) (p, typ ( P , t) + dt\- P , t)dj + > (-p, t) 



c ,[ +) ( P , t)cj-> ( P , t) - d { r) (_p, iK ( +> (-p, t) 
c , (+) ( P , t) c H ( P , t) - dj-> (-p, iK (+) (-p, t) 
4 +) (p, (p, *) - (-p, (-p, «) 



(23) 
(24) 
(25) 



dp (- Pl t)4 +) (-p, f) - 4 +J (p, <)c 2 - J (p, t) 



C+)/ 



„(-), 



(26) 



This representation can be named the complete quasi-particle one. The equations of motion for these amplitudes 
follow from Eqs. ifTBl) . (EH) 



cf > (P, t) = =■ A(p, *K +J ("P, *) ± MP, *)cr ; (P, t) + fl|f ; (p, t)cf J (P, f), 



At) 



7(±) 



-p,t) = ^A(p,i) C , (T) (p,i)± l u;(p,i)^(-p,i)+.9 l V (p,i)dP(-p,t). 



j(±) 



*(±) 



if)' 



(27) 
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The the matrix gij is defined as 



7(±) = 



±zeie2 






=Fieie 2 



£i(± 



e.i 



rn 
-A 



(28) 



where = (ei ± ie 2 )/-\/2. 

The transition to this representation from the initial one (a, b) is defined by the combination of the transformations 
© and Pljl. 



c (±) (p,t) = (7( ± )(p,t)a( ± )(p,t), 
d (±) (-p,t) = (7 (± '*( P ,i)b( ± )(-p,t), 



with non-unitary operator 



?7 (±) (p,i) = i? (±) -^(p.i) 



At) (t) (t) 

e l e 2 e 3 



.77631 — e 3 2 — e 3 3. 



(29) 



(30) 



The quantization problem is to be solved while taking into account the equation of motion l|27(l . It leads to the 
following non-canonical commutation relations 



[ct\p,t),c^\p',t)} = [d^(p,t),dl +) (p',t)]=Q^(p,t)Q^ } (p,t)S ppl , 
where the matrices Q^\p,t) are defined by the equations 



,(±) 



with the initial conditions 



lim Q[f\p,t) 



dij, 



(31) 



(32) 



(33) 



i.e. the commutation relations (|31|) transform to the canonical form only in the asymptotic limit t — > —00. The 
commutation relations l|31|) provide the definition of positive energy quasiparticle excitations with some time dependent 
energy reservoir of the vacuum. 



III. KINETIC EQUATION 



The standard procedure to derive the KE 3] is based on the Heisenberg-type equations of motion or l|27|) . Let us 
introduce the one-particle correlation functions of vector particles and antiparticles in the initial (a, 6)-representation 



F^(p,t) = (0 in \a(+\p,t)ai-\p,t)\0 in ), 
fUM = (0in|6L _) (-P.*)^ +) (-P.«)|0<n), 



(34) 



where the averaging procedure is performed over the in- vacuum state 2] . Differentiating the first one with respect to 
time, we obtain 



iV(p, t) = \ A(p, t) {f« (p, t) + Ftf (p, *)} , 
where the auxiliary correlation functions are introduced as 

F<±>(p,t) - (0 m |aW(±p,t)6( ± )(±p,i)|0 m ). 



(35) 



(36) 



The equations of motion for these functions can be obtained by analogy with Eq. (|35|l . We write them out in the 
integral form 



F^Hp,t) = ^ I dt'A(p,t') F^(p,t') + F^(p,t') 



1 



(37) 
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where 

0(p;Mo) - / dt'uip,?). (38) 



In Eq. l|37jl. the asymptotic condition F^\p, — oo) = (the absence of quasi-particles in the initial time) has been 
introduced. The substitution of Eq. i|T7|l into Eq. l(3l)|) leads to the resulting KE 



F^(p,t) = ^A(p,t) / dt'A(p,t')[F^(p,t') + F^(p,t')}cos[2 6(p-,t,t% (39) 

— oo 

This KE is an almost natural generalization of the corresponding KE for scalar particles ^| • 

Thus, the OR turns out to be an effective method for the diagonalization of the Hamiltonian in the Fock space. It 
is sufficient for the derivation of the KE Ij39(l . However, on this stage, there is a number of problems that are specific 
for the vector field theory: the energy is not positively defined, the spin operator has a nondiagonal form in the space 
of spin states etc., see above. This circumstance hampers the physical interpretation of the distribution function 
(|34|l . In order to overcome this difficulty, it is necessary to pass on to the complete QPR in which the system has 
well-defined values of energy, spin etc. The simplest way of derivation of the KE is the QPR, based on the application 
of the transformations l|29(l directly onto the KE (|39(l . 

A. Kinetic equation in QPR 

By analogy with the definitions l|34|l . let us introduce the correlation functions of vector particles and antiparticles 
in the complete QPR 

fik(p,t) = (0<„|^ +) (p l t)4" ) (P.*)|0<n), 

/<*(p,*) = (0i„\4-\-p,t)d { +\-p,t)\0 in ). (40) 

They are connected with the primordial correlation functions (134(1 by relations of the type 

f ik (p,t) = U+(p,t)U-Jp,t)F nm (p,t), (41) 

where F nm (p, t) is the "spatial" part of the tensor function F^(p, t) l|34|) (to, n — 1, 2, 3). 

To obtain the resulting KE in the complete QPR we differentiate the function fik(p,t) (|41|l with respect to time, 
and take into account the KE (|39H 



+ ^+\t)ui;\t)A(t) I *'A(^)C/ J ( + ) - 1 (^)^ K1 (^)[/ m „(t') + / m n(i')]cos20(p;t,i'). (42) 



i n 
t 

In comparison with non-Markovian effects of vacuum tunneling of scalar particles 0] , the considered case has its own 
characteristics related to the dynamics of spin twist. The Markovian limit (f mrl (t') — > fmn(t) in the integral part of 
r.h.s. of l)42[ll is admitted for rather slow processes and results in a significant simplification of the KE l)42[l. 

The system of the integro-differential Eqs. (|42|) can be reduced to a system of 27 coupled ordinary differential 
equations that is convenient for numerical calculations. We will not analyze here this rather complicated case and 
will restrict ourselves below to the consideration of a simple particular case having cosmological motivation. 



B. Deformation of the energy gap 



Let us consider the vacuum creation of vector bosons in the case when it is caused by an arbitrary time dependent 
deformation of energy gap, i.e. m = m(t) and A k (t) — 0. This is an isotropic case with Pf. = p^. and = 0. These 
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conditions are rather characteristic for different inflationary models of the preheating process (see, e.g., [l!j and the 
references given there). As a result, the KE l|42l) takes the following form 



fik{p,t) = - A m (p,t)[<W 3 fc(p,*) +<Wi 3 (p,t)] + 



1 



A(p,t) / ^'A(p,i')M^ 7 (p,t,^)[/ J 7(p,i') + / J 7(P^')]cos20(p;i,i'), (43) 



where 



and 5„ J 



uj{t') m{t) 
Lu(t) m(t') 



(44) 



As to be expected, the distribution functions f a p(p,t) and F a p(p,t) satisfy the same KE (|39|) for a = 1,2. The 
feature of the complete QPR becomes apparent only in the component of tensor distribution function f%k(p,t) that 
contains the preferred values of spin index i, k = 3. Let us select the KE for the diagonal components of the 
correlation functions (|40|l having a direct physical meaning as the distribution functions of the transversal (i — 1,2) 
and longitudinal components 



fi(p,t) = ^A(p,t) J dt'A(p,t , )[l + 2f i (p,t , ))co S 26(p;t,t'), 



(45) 



f 3 {p,t) = -2A m (p,i)/ 3 (p,i) 

t 

1 . , .,m 2 (t) 



+ o A (P>*)- 



2 (*) 



J dt'A(p, 0^5^ [ 2 /3(P> + Q(P> 0] cos 20(p; t, 0- 



(46) 



Here the shorthand notation /,j = /j has been introduced for the diagonal components of the matrix correlation 
functions 1)40(1. and we have 

A =—a-, A m = -A-^. (47) 

It is possible to show, that the distribution functions of the longitudinal (i = 3) and transversal (i = 1,2) components 
are connected by the relation 



/3(p,*) = Q(p,*)/i(p,i), 



(48) 



where Q(p,t) is the function occurring in the commutator of the creation and annihilation operators for the longitu- 
dinal bosons, 



[4 \p,t),ci + \p',t)} = Q(p,t)V, 



Q(p,t) 



exp 



-2 f A m {t')dt' 



m(t) u)(t ) 



1 2 



m(t ) u(t) 



(49) 



where m(to) and a; (to) are some initial values of m(t) and ui(i). 

Owing to Eq. (|48|1 . it is sufficient to solve the one equation l|45|) . We use now the well-known procedure of the 
reduction of the KE from the integro-differential form to the corresponding system of ordinary differential equations 
[20| in order to study the KE (|46|l numerically and to investigate the asymptotic behavior of its solutions for large 
momenta below in Sect. 



fk = -Aufc, iik = A(l + 2/fe) - 2uv k , 



v k = 2uju k . 



(50) 



Here u k and v k (k = 1, 2) are some auxiliary functions responsible for the different effects of vacuum polarization (see, 
e.g., |23). It can be shown by analogy with the scalar field case [T^l that this system is conservative and has the first 
integral of motion 



(1 + 2/ fc ) 2 - u\ - v\ 



(51) 
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if the particles are absent in the initial time, 

/fc(-oo) = Uk(-oo) = ffc(-oo) = 0. (52) 
The general initial condition for all diagonal components of the distribution function 



leads to the following requirement 



lim fk(t) = lim Uk(t) = lira Ufc(i) = (53) 

t — > — oo t — > — oc t — > — oo 



lim m(t) = mo, or lim m(t) = 0. (54) 

t — > — oo f — > — oo 



The main characteristic of the vacuum creation process is the total number density of vector bosons, 

2 oo oo 

n tot (t) = 2j2m(t) = ^ J p 2 d P [2h{p,t) + h{p,t)\=^ j p 2 dpfx{p,t)[2 + Q(p,t)l (55) 
1=1 o o 

where isotropy of the system was taken into account, p = |p|. The factor 2 corresponds to equal numbers of particles 
and anti-particles. As it will be shown in Sect. the integral (|55|l is convergent. 

Up to now, the time dependence of particle mass was brought in on the phenomenological level, without an indication 
of any concrete mechanism of its origin. As it was shown above, this is sufficient for the formal construction of kinetic 
theory of vacuum particle creation in the case of rather fast mass change. We will consider below the case when the 
time dependence of vector boson mass is defined by the conformal evolution of the universe. 



IV. VECTOR BOSON PRODUCTION IN THE EARLY UNIVERSE 



The description of the vacuum creation of particles in the time dependent gravitational fields of cosmological models 
goes back to Refs. |32> H3, 13 and has been recently reviewed, e.g., in the monographs [35|,|3(|. The specifics of 
our work consits in the consideration of vacuum generation of vector bosons in the conditions of the early Universe 
in the framework of a conformal-invariant cosmological model |14| , thus assuming, that the space-time is conformally 
fiat and that the expansion of the Universe in the Einstein frame (with metric g^) with constant masses m can be 
replaced by the change of masses in the Jordan frame (with metric g^ v ) due to the evolution of the cosmological 
(scalar) dilaton background field H,H3. This mass change is defined by the conformal factor VL(x) of the conformal 
transformation 

g^x) = n 2 {x)g^ v . (56) 
As mass terms generally violatate conformal invariance a space-time dependent mass term 

has been introduced which formally keeps the conformal invariance of the theory |2l| . In the important particular 
case of the isotropic FRW space-time, the conformal factor is equal to the scale factor, fl(x) — a(t), and hence 
mit) = a(t)m i, si where t is "the Einstein time" and m i, s is the observable present-day mass. Such a dependence was 
used, e.g., in Ref. [2| for the Robertson- Walker metric. On the other hand, the scaling factor a(t) is defined by the 
cosmic equation of state (EoS). For a barotropic fluid, this EoS has the form 

P P h = (7 - = & P h> (58) 

where p p h and e p h are phenomenological pressure and energy density (in the distinction from "dynamical" p and e 
(see Sect. 0) below), 7 is the barotropic parameter, c s is the sound velocity. The solution of the Friedman equation 
for such EoS leads to the following scaling factor 

o(t) - t 2/37 . (59) 

The kinetics of vacuum creation of massive vector bosons (Sects. II and III) was constructed in the flat Jordan frame 
with the proper conformal time t, which is necessary to introduce now in the Eq. I|59(l . The transition to the conformal 
time is defined by relation dt = dt/a(t). From this relation and Eq. (|59Jl follows 



t 



37 



37/(37-2) 

(60) 
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The substitution of this relation in the Eq. (|59|l leads to the mass evolution law in the terms of the conformal time 

2 

m(t) = (t/t H ) a m w , a = -, (61) 

07 — 2 

where tn = [(1 + a)!!]^ 1 is the scaling factor (the age of the Universe), H is the Hubble constant and mw — 80GeV 
is the W-boson mass. Let us write here the values of the parameter a for some popular EOS: 7 = 2, a = 1/2 (stiff 
fluids); 7 = 4/3, a = 1 (radiation); 7 = 1, a = 2 (dust). 



0,30 r 




1,0x10"" 3,0x10" 5,0x10" 1 2 

Time[s] P /m o 



FIG. 1: Time evolution of the particle number density with initial condition mo - to = 1 f° r a = 1/2 (left) and the corresponding 
momentum distribution at the time t 3> to (right). 
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FIG. 2: The dependence of the final particle number density on the initial time: a = 1/2 (left), a = 1/3 (right). 



Due to back reactions and dynamical mass generation during the cosmic evolution the detailed mass history remains 
to be worked out. The central question, however, is whether the number density of produced W-bosons could be of 
the same order as that of the CMB photons, ticmb ~ 465 cm~ 3 . If this question could be answered positively, the 
vacuum pair creation of W-bosons from a time-dependent scalar field (mass term) could be suggested as a mechanism 
for the generation of matter and radiation in the early Universe. The non-Abelian nature of the W-bosons could even 
imply consequences for the generation of the baryon (and lepton) asymmetry due to topological effects [8j . 

The numerical analysis of the Eqs. I|5l)|l is performed by a standard Runge-Kutta method on a one-dimensional 
momentum grid. As one can see from Fig. 1, the creation process ends very quickly and the particle density saturates 
at some final value. The momentum distribution of particles is formed also very early when m(t) « uiq and frozen in 
such form so that later on for times t 3> to most of the particles have very small momentum p -C m{t). The spectrum 
of created bosons is essentially non-equilibrium, hence we should continue further the analysis of relevant dissipative 
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mechanisms and other observable manifestations of the non-equilibrium state (e.g., CMB photons; in this connection, 
see, for example, p^l. 

The other interesting feature of the mechanism of particle production considered above is the qualitatively different 
form of momentum distribution when compared to the creation in the electric field. As is known, the electric field 
generates the pairs with monotonic momentum distribution and/(p = 0) = max in the fermion case and with 
non-monotonic distribution and f(p = 0) = in the boson case [2j. The creation mechanism connected with time- 
dependent effective mass modifies this situation to the contrary: the boson have a monotonic distribution (Fig. ^) 
and the fermions a non-monotonic one. The dependence of the corresponding final value of density on the initial time 
is shown in Fig. [5] The final density n\ of particles with spin projection ±1 reaches a maximum, when we let the 
initial time go to very early times, close to the birth of the Universe. However, in the same limit, the density / 3 of 
particles with spin projection zero grows beyond all bounds. The choice of the EoS changes drastically the quantity of 
the created particles, thus giving values which are too small (a = 1/2) or too large (a — 1/3) in comparison with the 
observed CMB photon densities. In order to improve this model, we should use an improved EoS, assuming that the 
barotropic parameter 7 charcterizing the evolution of the particle masses can change during the time evolution. Such 
a time-dependence could be induced by the back-reaction of the created particles on the scalar field. Furthermore, 
we could use another space-time model, e.g., the Kasner space-time |23j instead of the conformally flat de Sitter one. 
The main achievement relative to the earlier work @ is that in the™esent approach, there is no divergence in the 
distribution function, thus we do not need to introduce some ambiguous regularization procedure. 



V. EOS FOR THE ISOTROPIC CASE 



The relations for the energy density and pressure can be obtained from the energy-momentum tensor corresponding 
to the Lagrangian J2J 

Tn V = -d^v-advii - d v u* a d^u a - g^C. (62) 

Thus, we have e = (0 in |T O o|0 ir j) and p = l/3(0 in |T ii |0 i „). Let us substitute the decompositions JSJ and use the 
relations, which are a consequence of the spatial homogeneity of the system (in the discrete momentum representation) 

(0 in \a(+\p,t)al-\p',t)\0 in ) = (2^) 3 <5p, p ,iV(p,i), 
(0 in \b^(-p,t)b^(-p',t)\0 in ) = (2^) 3 <5p,p^(p,i), (63) 
(0 m | a ( ± )(p,t)fo( ± )(-p',t)|0 m ) = (2^) 3 ,5p,p,F(±)(p,£). 

As a result we obtain the following expressions for the energy density and the pressure (after the transition to the 
thermodynamic limit, dT = (2ir)~ 3 d 3 p) 

e(t) = -Jdru(p,t)[F£(p,t) + F£(p,t)], (64) 

+ (2p 2 + 3m 2 )[^(p,i)+F^(p,t)]}. (65) 

We perform now the series of the consecutive transformations of the functions Fff(p, t), Fjf(p, t) and F^ ± ^(p, t): the 
exclusion of the /i = component (according to Eq. and the transition to the complete QPR with the help of 

Eqs. I|30|) . Thus, we arrive at 



*£(P,*) = -i>(P'*)' ^(P.*) = -E/»(P.*). (66) 
1=1 1=1 

F^(p, t) = - £ jf \p, t) - V t), (67) 
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where fi(p,t), fi(p,t) and f^(p,t) are the corresponding functions in the complete QPR. The substitution of these 
relations into Eqs. I|64|) and Ht)5fl leads to the following expressions for the energy density and pressure 



3 

e(t) = 2]T / dTu>(p, t)ft(p,t), 
1=1 

= lltf : ^ 7 y{2p 2 / i (p,<)-[2^ 2 ( P ,t)+™ 2 ] 



(68) 



2p 2 
3m 2 



Taking into account the isotropy of the system, we obtain the EOS for the massive vector boson gas 



e(t) 
p(t) 



dT 



2 J drw(2 + Q)/i, 
[2p 2 (2 + Q)A] +S Pvac {t), 



where Sp vac (t) is the contribution in the pressure induced by the vacuum polarization, 



Sp V ac{t) = -- 



dT 



2uT 



1 



Ul. 



(69) 

(70) 
(71) 

(72) 



In order to prove the convergence of the integrals l|7U|l - (|72|l we investigate the asymptotic behavior of the solution 
to the system of Eqs. (I50|l . This system can be solved exactly in the asymptotic limit p ^> m for the case a — 1/2 
when it gets the form 



/ = 



in 



H 



1 



u 

V 



tt H p 2 

-±1.(1 

2pu 



2f) - 2pv, 



(73) 



The solution of i|73[l with the initial conditions (|53|l is 



f(p,t) 
u(p, t) 



sin 2 pit — t ) 
16(p/mo) 6 ' 
1 sin2p(t-t ) 
2 VM = 4(p/ mo )3 - 



(74) 



where niQ = m(to). The numerical investigation of the general Eqs. I|5(J|I shows that the basic features of the solutions 
(|74|l for a = 1/2 are conserved also for other a > 0. 

According to (|74|l . the particle and energy densities (Eqs. (|55|l and J7UJ), respectively) are convergent, but the 
vacuum polarization contribution to the pressure l|72() is divergent. Moreover, irrelevant fast vacuum oscillations of 
the pressure are observed here. Let us remark that such a behavior of the pressure for a plasma created from the 
is not a special feature of the present theory but is characteristic also for the models where an electron - positron 
plasma is created in strong, time-dependent electrodynamic fields as investigated in [24| . The standard regularization 
procedure of similar integrals with some unknown functions satisfying ordinary differential equations is based on 
the investigation of asymptotic decompositions of these functions in power series of the inverse momentum, l/|p|" 
(the procedure of n- waves regularizations [25|l. In the considered case, such a procedure is not effective because the 
solutions l|73|) has fastly oscillating factors ("zitterbewegung" ), the asymptotic decompositions of which lead to the 
secular terms. Therefore we regularize the pressure by a momentum cut-off at p = lOmg and separate its stable part 
by the time averaging 



< p >= 



1 



(t - 1 ) 



p{t)dt. 



(75) 



12 




FIG. 3: Time dependence of energy density and pressure at t 3> to with initial condition rm ■ ti = 1 for a = 1/2. 



The such " coarse graining" procedure was proposed in 26] in order to exclude the " zitterbewegung" from the descrip- 
tion of vacuum particle creation. In reality, the smoothing of these fast oscillations occurs due to dissipative processes 
that are not taken into account here. 

Fig. shows that the mean pressure remains negative and its magnitude becomes negligible in comparison with 
the energy density. At large times the energy density grows but the pressure stays very small, p ~ 0. The energy 
growth with the condition p ~ leads the conclusion mentioned in Sect. II VI that the massive vector boson-anti-boson 
gas created from the vacuum is a cold one (see also Fig. 1). It can be seen directly from Eq. JTHJ), that at large 
times e(t) ~ m(t)n(t), because of u>{t) ~ m(t). Let us remark also that such an EoS of the massive vector boson 
gas (e ^ and p ~ 0) corresponds to dust-like matter j23, which would characterize the evolution of the Universe 
during those stages when the vector boson gas is the dominant component of its matter/energy content. On the 
qualitative level, this conclusion is valid independent of the concrete choice of an EoS and, in particular, in the case of 
the dust-like EoS. In a subsequent work, we plan to obtain a formula of the type (|fil|) as a result of the solution of the 
Friedman equation with the EoS (|7UI) - (|72[) (such a procedure represents the back reaction problem) and to investigate 
self-consistently the production of vector bosons in the Universe. 



VI. SUMMARY 



The present work was devoted to the kinetic description of vacuum creation of massive vector bosons caused either 
by a time dependence of the mass or by the action of a non-stationary electric field. The statement of the problem 
is stimulated by modern cosmological problems related to the need for an explanation of the nature of the recently 
uncovered accelerating expansion of the Universe. The resulting KE (|42|l of non-Markovian type is obtained within 
a powerful non-perturbative framework and using the OR which provides an effective approach to the QPR, in the 
language of which the kinetic theory is constructed. We apply then this KE for the analysis of the important particular 
case of an isotropic gas of vector bosons with a time dependent mass which can be justified on the basis of a conformal- 
invariant scalar-tensor theory of gravitation. We show that the kinetic theory leads to a reasonable density of vector 
bosons in an early period of the Universe evolution which is sufficient to explanation the present density of CMB 
photons. 

The obtained results constitute a foundation for the subsequent investigation of the dynamics of vector bosons 
created from the vacuum (the equation of state, the long wave-length acoustic excitations, the back-reaction problem 
etc.). It is necessary to underline that we have considered here the single mechanism of mass change induced by 
the conformal expansion of the Universe. It was necessary to switch on the mass mo (to) at some arbitrary initial 
time to- Thus, for the construction of a more consistent theory, one should eventually take into account the inflation 
mechanism of mass generation acting during an earlier period the Universe evolution |l9t l28l l29l l3("| . 

It is very interesting to investigate the generation of elementary particles of different masses. Eq. (|61|l is valid for 
all elementary particles independent of their inner symmetry class. Then the chemical composition of created matter 
and the EoS of the Universe must be fixed rather precisely and can be subject to experimental verification. 
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